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ABSTRACT 


This  paper  presents  a characterization  of  the  solutions  of  a special 
quadratic  program.  This  characterization  is  then  used  in  the  development 
of  an  efficient  algorithm  for  this  class  of  problems. 
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II.  CHARACTERIZATION  OF  SOLUTIONS 


Let  us  rewrite  (1)  - (3) 

in  a slightly  different  form  as 

follows: 

min  J x"Dx  - a^x 

(A) 

s.t.  ZjXj  - c - 0 

(X) 

(5) 

xs  ' b s - 0 

(V 

(6) 

- Xj  < 0 

<Vj)’ 

(7) 

where  X,  , and  are  the  Kuhn-Tucker  multipliers  associated  with  the 
three  types  of  constraints.  The  Kuhn-Tucker  conditions  for  (4)  - (7)  may 
be  stated  as  follows: 


djXj  - a^  + Uj  - Vj  + X ■ 0,  (for  all  j) 

(8) 

uj(Xj  - bj)  “0,  (for  all  j) 

(9) 

VjXj  * °’  (f°r  a11  ^ 

(10) 

Uj  ,Vj  _>  0,  (for  all  j) 

(11) 

plus  (5),  (6),  and  (7), 

where  dj  is  the  jth  diagonal  element  of  D.  Consider  the  following  solution 
as  a function  of  X. 


For  any  selection  of  A,  the  above,  solution  clearly  satisfies  (6),  (7),  and 
(11).  We  now  show  that  this  solution  will  also  satisfy  (8),  (9),  and  (10). 


Proposition  1. 


The  solution  given  by  (12)  satisfies  (9). 


Proof.  Case  1:  a.  - A - d,b.  < 0 »>  u.  « 0 *>  u.(x.  -b.)  =0 
j j j j i i i 

a,  - A 


Case  2:  a - A - d b > 0 =>  > b.  =*>  x.  = b.  =*>  u.  (x  - b.)  = 0. 

i i 1 - dj-j  3 3 j'j  j" 


This  completes  the  proof  of  proposition  1. 


Proposition  2. 


The  solution  given  by  (12)  satisfies  (10). 


~ A 

Proof.  Case  1:  a.  - A < 0 ■>  — < 0 »>  x = 0 =>  v x = 0. 

J d.  j j J 


Cas?_2,  aj  - X > 0 ■>  vj  - 0 ->  vjxJ  - 0. 


This  completes  the  proof  of  proposition  2. 


Proposition  3. 


The  solution  given  by  (12)  satisfies  (8). 


Proof:  Case  1: 


a - A 


*T- bJ ->  ' V aJ ' x ' Vj  - and  ai ' x * °- 


*j  ' x ‘ Vj  i o ->  Uj  - - x - djbr 

3j  . X > 0 ->  Vj  . 0.  Thus  Yj  - + «J  - VJ  + x 

" Vj  ‘ aJ  + *J  ' * ' Vi  ‘ 0 + J ' °' 


a.  - X a.  - X 

o < J__  < bj  .>  X}  . -L— , 3j  - X - d.bj  < 0,  and  a.  - X > 0. 

aj  - x - djbj  < 0 •>  Uj  . 0. 

*J  - A > o ->  v3  - 0.  Thud  djXj  - + -j  - vj  + x 


■■(V, 


-aj+0-0  + X = 0. 


Case  3: 


a i - X 

— J- < 0 =>  x.  = 0 and  a.  - X < 0. 

di  ~ 3 3" 

Sj  -x<0->aj  -X-djbj  <0  andVj-X-aj. 

bj  - X - djbj  < 0 ->  u.  - 0.  Thus  djXj  - a.  + Uj 


VJ  + X 


■ 0 - aj  + 0 - (X  - Sj)  + x - 0. 

This  completes  the  proof  of  proposition  3. 

Hence  to  solve  (1)  - (3)  one  need  only  find  the  appropriate  X such  that  (5) 
is  satisfied. 


"Hf 

Then  we  must  find  X such  that  g(X  ) * c.  Note  that  since  d^  > 0 and  > 0, 
x^(X)  may  be  expressed  as  follows: 


Clearly  each  (A)  is  piece-wise  linear  and  monotone  nonincreasing.  Since 
the  sum  of  such  functions  preserves  this  property,  g(A)  is  piece-wise  linear 
and  monotone  nonincreasing.  A typical  g is  illustrated  in  Figure  1. 


Figure  1 About  Here 
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III.  ALGORITHM 


Suppose  there  are  n (X) , then  the  breakpoints  for  the  piece-wise 


linear  function  g(A)  occur  at  the  2n  points  a^  - d^b..  and  a^  for  j = l,  . ..,n. 


Let  y, , . . . , y„  denote  these  breakpoints  where  y.  < y„  < ...  < y0  . Then 
1 zn  l — z — — Zn 


for  X £ y^,  g(A)  = Z^b^  an^  ^or  ^ — ^2n’  = Hence  (1)  “ (3)  has  a 


feasible  solution  if  and  only  if  0 < c < I.b.. 

— — J J 


We  now  present  an  algorithm  for  obtaining  the  value  X , such  that 
* * 


g(A  ) = c.  The  procedure  consists  of  a binary  search  to  bracket  X 
between  two  breakpoints  followed  by  a linear  interpolation. 


ALGORITHM  FOR  A* 


7 nttializalion 


If  c > Ejbj  or  c < 0,  terminate  with  no  feasible  solution;  otherwise,  set 


SL  1,  r ■*-  2n, 


L - Zjbj  and  R - 0. 


Teit  Foa.  Bracketing 


If  r - H = 1,  go  to  4;  otherwise,  set  m *■  [(£  + r)/2]^,  where  [k]j  is 


the  greatest  integer  <_  k. 


2.  Compute  New  Value 

Set  C MAX  |m1N 


~ ym  \ 

rV-  v ■ °. 


Update 


X 

If  C » c,  terminate  with  X «- 


m 


If  C > c,  set  Jl  m,  L C,  and  go  to  l. 
If  C < c,  set  r «-  m,  R ■*-  C,  and  go  to  l. 
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